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Abstract. Supermassive disks are thought to be precursors of supermassive black 
holes that are believed to power quasars and exist at centers of galaxies. Formation 
scenarios of such disks are reviewed and it is argued that gas dynamical schemes are 
favourable compared to stellar dynamical schemes which could however be important 
feeding mechanisms for the growth of the black hole. A new self-similar model of a 
collapse of a self-gravitating disk due to radiation induced stresses applicable to two 
different situations of radiative viscosity and Compton drag is presented. The collapse 
timescale purely due to radiative viscosity is found to be a fraction of Hubble time, 
Ty ~ (Ttc I {mpG){Ledd I L) ^ 6 x lO'^yrs is slow and probably magnetic fields play an 
important role before general relativistic effects take over. A model of self-gravitating 
disk collapsing due to Compton drag by the Cosmic Microwave Background is also 
presented which is found to be effective at redshifts 1400 > z; 300. It is proposed 
that the small < IO^'Mq objects that form by this mechanism by z ~ 20 can merge and 
coalesce by dynamical friction to form the high redshift quasars seen. Supermassive 
stars which are systems (and could be end products of a supermassive disk phase) en 
route to the final collapse are also briefly reviewed. 

Keywords : Black Holes -Formation, Supermassive Stars, Radiation Hydrodynamics, 
Cosmic Microwave Background. 



1. Introduction 

There seems to be increasing evidence that supermassive black holes are at the centers of galaxies. 
Dynamical searches indicate the existence of massive dark objects (MDOs) in eight systems and 
their masses range from 10^ - lO^^M© (Kormendy & Richstone 1995). Although this study does 
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not confirm that the central objects are supermassive black holes, it has been inferred that the 
central mass is contained within 10^ Schwarzchild radii. On an average, the black hole mass is 
a fraction, 10"^ - 10"-', of the total mass of the galaxy and of order 10"^-'' of the bulge mass 
(Wandel 1999). Recent observations show a strong correlation between the black hole mass, Itl/,, 
from stellar dynamical estimates, and the velocity dispersion of the host bulges (nVj oc cr"; where 
a is reported to be in the range 3.5-5; eg. Ferrarese & Merritt, 2000). 

The presence of quasars at high redshifts tells us that galaxy formation had proceeded far 
enough for supermassive black holes to form in the standard picture (Rees 1984). A detailed 
model of formation of these objects should address the issues of supernovae feedback from star 
formation and the mechanism of efficient angular momentum transport in order to explain the 
massive active nuclei as early as z = 6. In the case of MDOs, there is a need to explain the 
compact sizes of 10 - 100 pc that are imphed from dynamical studies. 

2. Formation scenarios 

We first discuss the possibility that dense star clusters feed a seed black hole which grows initially 
by accreting the gas resulting from tidal disruption that occurs at the tidal radius, Rj ^ {6X&hl^P*)- 
Subsequently, when the Schwarzchild radius,/?^, exceeds this, the black hole then grows by 
swallowing low angular momentum stars whose pericenters he within Rs. The timescale for 
growth during the gas release through tidal disruption suffers from serious drawbacks- the ac- 
cretion rate is given by (see eg. Hills 1975), M = {An^GXAhRj I a)ps, assuming a Maxwellian 
distribution of velocities with dispersion cr. At this rate, the growth timescale goes as 1.7 x 
10"(p.J10'^Mopc"^)(m/,/Mo)'^^(cr/100km/sec)yrs which is large even for a dense cluster (the 
alternative of a seed mass of 1000 Mq for the black hole is unhkely from theories of stellar evo- 
lution). Such models then beg the question of how clusters of densities ~ lO^Mo/pc^ form in the 
first place. 

Now consider the situation where the growth proceeds by swallowing of stars. Swallowing 

rate of stars may be estimated in the following way. Stellar orbits diffuse by two-body relaxation 
toward lower angular momentum orbits until they enter a small loss-cone of semi-aperture 9c - 
{tdynltR)^^^ (Frank & Rees 1976, Lightman & Shapiro 1977) where is the relaxation time and 
tdyn is the dynamical time. The resulting swallowing rate is, M ^ Ncnit /(tRln(2/0c)) ~ mt/tdm 
which is not rapid enough in most cases. The key point is that in this scenario, one must assume 
a extremely dense and massive cluster. 

Recently, there have been proposals (Volonteri, Haardt and Madau 2003, Wyithe and Loeb 
2002) motivated by quasars discovered at z « 6 that these objects have are assembUes of smaller 
IOOMq objects that collapsed at z ~ 20 from high-cr density fluctuations. The former model 
invokes dynamical friction of merging systems to sink to the center as larger halo objects and 
involves Monte-Carlo computations based on halo merger trees from cosmological simulations in 
a A CDM cosmology. The model appears to yield the desired results of luminosity functions. This 
promising model for forming quasars at high red-shifts is worth exploring further to ascertain how 
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massive the seed black holes need to be to explain the high redshift quasars and the efficiency of 
mergers. Observations of ultra luminous X-ray sources (ULXs) by (Colbert & Mushotsky 1999, 
Makishima et al. 2000) in nearby galaxies seem to indicate that seed black holes of intermediate 
mass of a few hundred Mg are possible. 

There are good reasons to think that supermassive gaseous objects are remnants of a galaxy 
formation process. Mangalam(2001) presented a detailed physical model wherein protoquasars 
(or MDOs) form from a magnetized accretion of a collapsed disk, the properties of which are 
obtained taking into account supernovae feedback in a virialized halo. Significant star formation 
and supernovae activity occurs after the cloud, which is spun up by tidal torques, contracts to a 
radius where self-gravity is significant. The model is composed of the following stages- 

1. The formation of a gaseous disk with a radial extent of about a kpc, in a host galaxy as 
Umited by supernovae feed back. The range in halo mass for a given redshift that still 
retains the hot gas was calculated. 

2. In previous work, gravitational instabilities in the disk was considered as the main source 
of viscosity. Justification was made for a magnetic viscosity and the estimated accretion 
rate turns out to be significant. The collapse of the disk was calculated with a generalized 
viscosity prescription (which includes the individual cases of magnetic, a and self-gravity 
induced instabiUties, under a halo dominated gravitational potential into a compact central 
region at rapid rate of about a yr ' . A self-gravitating magnetized disk solution for this 
central object that collapses to a seed black hole in 10'' yrs, was calculated. 

3. The impUcations for quasar luminosity functions and the time delay between collapse and 
virialization is considered in Mangalam(2003) and is based on the mass limits from cooling 
considerations in Mangalam(2001). 

3. Collapse of supermassive disks 

Here we calculate the collapse of self-gravitating compact mass that takes into account radiative 
stresses, which is a Newtonian 1.5 dimensional version of a quasi- spherical relativistic collapse 
currently under investigation. A particular application can be made to disks collapsing under 
angular momentum transport by radiative drag due to CMBR at high redshift. Another application 
is to estimate collapse time scale due to radiative viscosity after sufficient accretion of mass into 
a compact region of radius ro, typically of the order of a hundred parsecs containing a mass 
of lO^M©. The problem of self-gravitating accretion flow is complicated by the coupling of 
Poisson's equation to the momentum and continuity equations. Clearly, its evolution has to be 
treated difl'erently from ffie case of a prescribed background potential. 

We consider a simpler model of disk where the self-consistent density distribution wiffi a 
gravitational potential that is entirely due to self-gravity, is of the Mestel form S(r, t) - v^(t) / {2nGr), 
where the time dependence appears only in the rotational velocity. Taking - VQxit) and 
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r = roxMx where Vq = GMdro, and is the mass out to vq. We see that by assuming a 
self-similar evolution of the disk, the mass out to a given x should be independent of t and hence 
it follows that;^i = x~'^ and ^ = ^cAf^/JC where = Mdilnr^). From the continuity equation, 

rdfi: + drifLvr) = (1) 

we find 

Vr = -^x'^xr^x. (2) 
Substituting this and the self-similar forms given above into the angular momentum equation, 

S5,v^ + ^(yrlr)dr{rv^) = il/r^)drir^Tlr^) (3) 

we obtain 

3 2 

= -^Vo^croX^X, (4) 

which is independent of x. So far no specific viscosity mechanism has been invoked- the form of 
Tlr<p above is necessitated by the prescription of a Mestel disk. If a stress due radiative viscosity 
is assumed, Ilr^ = r]yrH^, where the coeflicient of radiative viscosity (Misner 1968, Weinberg 
1971) rjy - (S/27)(ey)/io-Tnec), where Ej is the photon energy density, H is the half thickness, 
ctt is the Thomson cross-section and is the electron density. From the energy dissipation 
condition, the heat flux is given by 
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taking only the relevant component of the heat flux, 

., c 
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for Thomson scattering opacity. It follows that the half thickness, H oc xxifY^- Realistically the 
1.5-D assumption breaks down. Nevertheless, we can push our model to get some estimates. It 
follows that rjy oc Cylpc^ oc x. Since it is radiation dominated, we assume a polytrope of index 4/3 
and obtain 

\ryx'-X=x"\ (7) 
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which leads to the solution 



^(0 = (-.l) , (8) 

where ;(f(0) = 1 was taken as the initial condition. The factor eylipc^) can be estimated by 
calculating the luminosity due to the heat flux by taking Cy - 3py. The collapse timescale is then 
given by 

Ty- — '6xlOVs (9) 

mpU in L ro 

where = GMdc^, the gravitational radius; the fiducial value taken here corresponds to a 
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situation when the system is sufficiently compact and radiating at a tenth of eddington luminosity. 
Clearly, the model is not strictly valid when it is relativistic. The toy Newtonian model of self- 
similar, self-gravitating, collapsing due to radiative viscosity yields a collapse timescale, Ty, and 
suggests that the final phase of a Newtonian radiation dominated collapse to a black hole is 
slow. This example emphasizes the importance of other viscosity mechanisms like magnetic 
fields before destabiUzing GR effects take over; a detailed model is currently under study. 

The above model can be easily adapted to the case of coUapse due to radiation drag at high 
redshift; the corresponding angular momentum equation can be written as 

IStV^ + 'L{Vrlr)dr{rv^) = -kLv^, (10) 

where K{t) = ^ey{t)aT licMp) is the coefficient of Compton drag in a completely ionized plasma 
and Cyif) = aT^it) is the CMBR energy density. We obtain the solution = exp{J Kdt/3), 
and hence x\(t) - exp (-(2/3) j Kdt); further using z = -//()(! + z)^^'^ appropriate for a matter 
dominated era, the collapse factor due to Compton drag is given by 



l/Xiiz) = exp[i|//oi— cr^r^ [(1 + Zif - (1 + Zff^]] = exp 



'2K(zdt(zd 



(11) 

which shows that the e-folding time in the angular momentum at z, ^ 1400, is initially shorter 
than the Hubble time by two orders of mag nitude- Kizdtizi) = 2.2x10^ ((1 + Zi)/I400f'^, aresult 
that is consistent with Loeb(1994) who computed a spherical model of a cloud in an expanding 
background (here it is assumed that the disk has turned around and collapsed). Therefore CMBR 
drag is effective at redshifts 1400 > z> 300 and it is possible to collapse smaller mass clouds 
< IO^Mq. 



4. Supermassive stars 

Some of the proposed scenarios are envisioned to lead to a build up of a supermassive star at 
relativistically compact scales as an intermediate stage of the evolution before the gravitational 
instability sets in and a rapid final collapse to a supermassive black hole ensues. 

Supermassive stars (SMS) are equilibrium configurations that are dominated by radiation 
pressure (the luminosities are nearly at the Eddington limit) and can have masses between IO'^Mq 
and about lO^M©. They are expected to be fully convective (Loeb & Rasio (1994) give a for- 
mal argument that radial entropy has to develop eventually which drives convection), isentropic 
and their structure can be well described by a Newtonian polytrope with y = 4/3. The en- 
ergy per nucleon in an SMS is given by the radiation entropy in units of the Boltzmann con- 
stant (S/k ~ 0.94(M/Mo)'^^ where M is the mass of the star (see eg. Zeldovich & Novikov 
1967). The evolution proceeds on a Kelvin-Helmholtz timescale and is driven by loss of energy 
and in the case of rotating SMSs, loss of angular momentum through mass-shedding. Pres- 
sure contributions from plasma components raise the adiabatic index of the equation of state. 
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r = 4/3 + jS/6 marginally above the critical value 4/3 where yS is ratio of gas to radiation pres- 
sure. General relativity leads to the existence of a maximum for the equilibrium mass (as a 
function of the density) and gravitation instability sets in when F falls below a critical value 
r < ~ (2/3)(2 - 577)/(l - lif) + l.l2Rs/R where 77 = T/\W\ is the ratio of the rotational 
to the gravitational potential energy and R is the radius of the star (Misner, Thome & Wheeler 
1973). If the plasma contribution is not enough, then the star shrinks during the evolution; ro- 
tation can however hold up the collapse if 77 is above a critical value. The SMS evolves on a 
Kelvin-Hehnholtz timescale tKH = \Ec\/Ledd - 10^/(M/Mo)yrs where Ec = 3 X lO^'^ergs is the 
equilibrium energy at the onset of gravitational instability. Typical SMSs with M ~ IO^Mq have 
a Ufetime of a 1000 yrs. Rotation can appreciably stretch their equihbrium evolution; Baum- 
garte and Shapiro (1999) found a lifetime independent of stellar mass of (9 x lO'^s) and that key 
non-dimensional ratios, R/Rs, t], and Jc/(GM^) for maximally and rigidly rotating polytropes are 
independent of the mass, spin or radius of the star. Driven by radiation and angular momentum 
loss through a mass-shedding sequence, the SMSs collapse leads to a formation of a black hole 
through an explosion powered by hydrogen burning in CNO cycle and associated with gigantic 
release of neutrinos. However, this neutrino flux and resulting background from such sources is 
weak for even the new generation neutrino detectors like Super Kamiokande; but the possibility 
that about 10% of the baryons are locked in SMSs at z < 1 (Shi & Fuller 1998) can potentially be 
ruled out. 



5. Summary and discussion of observational discriminants 



Based on the arguments here models of monolithic formation and collapse (Mangalam 2001, 

2003) of supermassive disks can explain the mass and redshift distribution of black holes of z< 6. 
However plausible merger scenarios (with dynamical friction) of smaller < lOOMg objects that 
form around z - 20, probably through the direct collapse and supermassive star route are needed 
for highest redshift quasars. In the case of self-similar self-gravitating contraction due to radiation 
viscosity, the timescale turned out to be a fraction of Hubble time, Tj ~ o-Tc/inipG) = 6 x lO'^yrs; 
clearly, one needs to take into account GR effects which is a future goal. The case of Compton 
drag yielded a timescale that is two orders of magnitude smaller than the Hubble time and is 
effective in the range 1400 > z> 300. Hence, it is possible to collapse smaller mass clouds 
< 10^ Mq. It worth investigating the statistics of such collapsed objects and whether mergers can 
produce the high redshift quasars in the picture above. The luminosity function needs to be 
more precisely determined to help distinguish between the models and the resulting stellar cusp 
profiles of the merged systems are likely to be different. The observed relation, VHf, oc a^, with 
a = 4 - 5 can be explained on the basis of the following physical arguments of saturation of 
black hole mass in (Silk & Rees 1988, Wyithe & Loeb 2002): the black hole mass saturates 
when luminosity impedes further accretion; ie, Ledd B.E./tdyn, where the gravitational binding 
energy scales as M^/R resulting in a = 5. Alternatively, though in a similar vein, the fraction 
of stars in an isothermal distribution that is captured by a Schwarszchild black hole is given by 
fir) ^ (7cap/(2crr))^, where /cap = 4GM/c is the maximum angular momentum for capture; this 
translates into a energy flux (evaluated near the radius of influence, r^) oc p(rh)r^crf(rh) which 
again results in a = 5 (Zhao, Haehnelt & Rees 2002). 
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In models of disk contraction that depend on self-gravity induced instabilities, the accretion 
is effective only upto the point when the Keplerian potential dominates over the gravity of the 
disk, which implies final black hole masses III/, > IQt'MqMcb where Mjq is the disk mass in units 
of 10^ Mq (Mangalam 2001, Loeb & Rasio 1994). This has direct impUcations for the seed mass. 
In conclusion, our understanding of the process is only begiiming and there several promising 
ideas that need to be further explored and more observations tests are required for discriminating 
amongst the models. 
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